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Bedford,  Massachusetts,  under  Contract  AF  19(628)- 
3852.  The  work  was  initiated  by  Mr,  Walter  Rotman 
and  monitored  by  him  and  Mr.  John  F.  Lennon. 

The  investigation  leading  to  the  results  re¬ 
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report  was  written  by  Mr.  Joseph  J.  Rossi  and  Mr. 
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ABSTRACT 


This  report  outlines  method  of  creating  a  win¬ 
dow  for  radio  communication  through  a  re-entry  plas¬ 
ma  sheath.  The  results  of  non- equilibrium  numerical 
calculations  illustrate  that  the  injection  of  a  coolant 
gas  will  reduce  the  plasma  sheath  temperature  but  will 
not  reduce  the  free  electrondensity  or  plasma  frequency. 
Plasma  electron  density  can  be  reduced,  however,  by 
attachment  to  electrophilic  gases.  The  chemical  rates 
of  attachment  are  computed;  and  criteria  are  developed 
which  specify  the  amount  of  electrophilic  gas  required 
to  reduce  the  electron  density  a  specified  amount  in  a 
given  time.  As  a  practical  application,  a  gas  injection 
system  for  the  total  alleviation  of  S-band  blackout 
throughout  the  re-entry  of  a  Trailblazer  II  vehicle  would 

have  a  total  charged  weight  of  S.  4  lbs. ,  a  volume  of  715 
3 

in,  ,  and  would  discharge  0.  Z8  lbs.  of  SF^  during  re¬ 
entry.  With  a  pulsed  injection  system,  an  even  smaller 
system  can  be  used. 
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1.  INTRODUCTION 


Plasma  blackout  of  radio  communications  during  re-entry  is  a 
problem  so  formidable  that  it  has  become  an  accepted  phenomenon  during 
the  Mercury  and  Gemini  programs.  Blackout  has,  so  far,  been  tolerated  on 
the  assumption  that  the  few  minutes  of  lost  communications  are  not  crucial. 

In  future  missions  such  as  Apollo,  however,  communications  losses  during 
re-entry  will  be  quite  critical  since  the  blackout  will  persist  for  tens  of 
minutes  and  could  therefore  threaten  the  safety  and  control  of  the  vehicle 
during  this  maneuver. 

Essentially,  radio  blackout  is  caused  by  the  presence  of  large 
concentrations  of  free  electrons  in  the  flow  field.  Figure  1  illustrates  the 
attenuation  that  occurs  around  a  blunt  body  from  the  flow  of  ionized  air;  the 
air  that  is  ionized  at  the  front  of  the  vehicle  expands  around  the  body  and 
cools.  If  the  air  could  at  all  times  be  in  thermodynamic  equilibrium,  the 
attenuation  would  be  reduced  as  illustrated  by  the  bottom  trace.  Unfortunately, 
however,  the  numerous  molecular,  atomic,  and  electronic  species  found  in 
air  react  together  at  a  finite  rate  which,  in  hypersonic  flight,  is  too  slow  to 
keep  up  with  the  rate  of  density  expansion  around  the  body.  Equilibrium 
cannot  thus  be  maintained,  and  attenuation  occurs  at  a  level  somewhere 
between  the  two  top  traces  of  frozen  and  non-equilibrium  flow.  The  extreme 
level  of  attenuation  shown  by  the  top  line  occurs  when  no  specie  change 
occurs  and  the  flow  is  considered  frozen. 

Creating  a  transmission  window  through  the  plasma  sheath  around  a 
re-entry  vehicle  requires  a  technique  for  reducing  the  electron  density  in  the 
regions  of  the  air  flow  surrounding  the  antenna.  The  obvious  way  to  reduce 
the  electron  density  is  to  cool  the  air  flow,  either  by  aerodynamic  shaping  or 
by  coolant  injection.  Unfortunately,  the  flow  can  be  cooled  without  an 
accompanying  reduction  in  the  plasma  electron  density  since  the  inherent 
electron  recombination  rate  is  too  slow  for  thermodynamic  equilibrium  to  be 
achieved  anywhere  near  the  point  of  coolant  injection.  To  speed  up  the 
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recombination  rate  an  additive  must  be  found  to  serve  as  a  new  "sink"  for 
electrons;  one  possibility  is  the  removal  of  electrons  by  electrophilic  gas 
attachment. 

The  best  known  electrophilic  gases,  the  halogens  —  chlorine,  fluorine, 
iodine  and  bromine  —  have  a  large  affinity  for  electrons.  In  fact,  equilibrium 
calculations  show  the  halogens  to  be  capable  of  reducing  the  plasma  electron 
concentrations  several  orders  of  magnitude  when  the  plasma  has  cooled  below 
3000°K,  Furthermore,  halogen  gases  are  commercially  available  in  many 
molecular  compounds  such  as  SF^,  CCl^,  Freon- 14,  Freon- 12  and  C^,  and 
these  compounds  are  themselves  electrophilic  and  can  easily  be  injected  into 
the  plasma,  some  under  their  own  vapor  pressure. 

Upon  injection  into  the  plasma,  electrophilic  gases  attach  electrons 
by  one  of  two  processes.  (1)  direct  three  body  attachment  between  the 
,  halogen  atom  and  electron,  or  (2)  dissociative  attachment  between  the 
halogen  molecule  and  electron.  Which  reaction  is  dominant  depends  on  the 
flight  altitude  and  location  of  the  injection  point  in  the  flow.  For  injection 
into  highly  expanded  flow,  the  frozen  electron  density  will  be  reduced  by 
dissociative  attachment,  and  injection  at  the  stagnation  point  will  favor  three 
body  direct  attachment. 
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2.  INJECTION  SYSTEM  CONCEPT 


To  cool  a  supersonic  flow  by  means  of  jet  injection  of  electrophilic 
gas,  the  maximum  lateral  penetration  must  be  achieved  with  as  little 
disturbance  to  the  main  flow  as  possible.  The  specific  method  adopted  is  to 
expand  the  injected  fluid  to  a  high  Mach  No.  and  to  a  pressure  which  balances 
the  pressure  of  the  external  stream.  This  correctly  expanded  jet  is  directed 
across  the  main  stream  so  that  the  momentum  carries  the  jet  fluid  well  out 
from  the  orifice  before  the  fluid  is  swept  downstream.  The  fluid  spreads 
into  the  stream  by  turbulent  mixing,  and  then  gradually  mixes  with  the  external 
stream  to  fill  the  shock  layer  adjacent  to  the  radio  antenna  as  shown  in 
Figure  2. 

Multiple  jets  are  used  so  that  the  airstream  is  not  blocked;  if  the 
airstream  were  blocked  a  strong  standing  shock  wave  would  be  set  up  in 
front  of  the  field  of  jets.  A  plan  view  of  the  jets  is  shown  in  Figure  3.  Each 
tiny  gas  jet  creates  a  shockwave  in  front  of  it  which  heats  up  only  a  small 
portion  of  the  oncoming  airflow,  most  of  which  passes  around  the  jet 
unchanged,  and  is  gradually  mixed  into  the  jet  flow  containing  the  electro¬ 
philic  additives.  In  this  way,  the  jets  do  not  alter  the  basic  airflow  while 
adding  the  electrophilic  gases. 

2.1  Jet  Size 


The  MITHRAS  technique  (Ref.  1)  for  designing  the  individual  jet  size 
and  Mach  No.  to  achieve  the  desired  penetration  is  presented  below.  The 
behavior  of  the  flow,  the  mixing  and  the  penetration  distances  will  depend  on 
injection  Mach  No.,  angle,  fluid  density,  and  free  stream  properties. 

In  Figure  4,  a  circular  jet  is  represented  by  an  expanding,  curving  tube 
in  the  external  stream.  The  tube  will  contain  more  and  more  external  stream 
air  entrained  by  turbulent  mixing  as  it  moves  outward.  Near  the  orifice,  the 
upstream  or  leading  edge  will  act  like  a  bluff  obstacle  in  the  main  flow, 
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and  there  will  be  a  normal  shock  in  the  flow  before  it  spreads  around  the  jet 
tube  and  closes  behind  it.  Actually,  part  cf  the  jet  fluid  will  be  torn  away 
from  the  tube  and  carried  downstream,  thus  distorting  the  circular  form  and 
reducing  the  total  mass  of  the  jet.  The  impact  of  the  free  stream  on  the  jet 
tube  will  cause  it  to  start  curving  downstream.  As  the  jet  progresses,  it 
picks  up  mass  by  entraining  external  stream  air  while  maintaining  a  constant 
momentum.  Since  the  jet  will  be  curving,  a  centrifugal  force  will  appear 
in  the  force  equations.  As  the  jet  tube  turns,  it  presents  a  different  angle 
to  the  external  stream,  thus  changing  the  impact  which  the  external  stream 
exerts  on  it. 

The  jet  penetration  is  determined  from  a  solution  of  the  force 
equations,  and  the  growth  rate  of  the  mass  flow  in  the  jet  is  given  by  Ricou 
and  Spalding  (Ref.  2).  Solutions  of  the  force  equations  are  given  by  Adams 
and  Rossi  (Ref.  1). 


y/d  =  a  J  sin  <J  (s,c)ds 

(2.1) 

x/d  =  a  f  cos  a  (s,~c)  d¥ 

(2.2) 

where 

M.  sin  a 

(2.3) 

- 2-47  (g) 

(2.4) 

7  -  (s/d)  (!) 

(2.5) 

Solutions  for  jet  penetration  (Figures  5,  6  and  7)  were  obtained  for 
the  injection  of  SF^  but,  unfortunately,  there  are  no  previous  experiments 
with  which  to  compare  the  present  findings  and  theory.  Several  reports  are 
available  on  penetration  from  underexpanded  sonic  jets,  and  References  3  and 
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4  present  height  to  diameter  ratios  of  between  10  and  50  for  incorrect 
expanded  flow.  As  already  mentioned,  though#  the  injection  of  electrophilic 
gases  requires  correctly  expanded  flows  so  that  the  jet  will  be  as  small  as 
possible  to  not  form  a  strong  shock.  Therefore,  wind  tunnel  experiments 
must  be  conducted  to  observe  the  penetration  of  correctly  expanded  jets. 


Jets  can  be  designed  for  a  specific  application  by  developing  design 
charts  from  the  numerical  results  presented  in  Figures  5,  6  and  7.  The 
jet  mass  flow  rate  of  SF^  through  a  sonic  or  supersonic  nozzle  can  be 
specified  as 


1 . 09  P.  A* 


(2.6) 


where  A  is  the  jet  throat  area  and  rhgp  is  the  jet  mass  flow  rate  of  SF^. 

6 

The  jet  exit  diameter,  dgx#  can  be  obtained  by  specifying  the  ratio, 

mQir  /P  . 

SF6  ex 


m 


SF/ 


ex 


1.09  a-  dex  /  A*  \  Pt 

"~r“  >Fo  \Aex‘  ^ 


(2.7) 


where  A  and  P  are  the  jet  exit  area  and  pressure.  The  ratios 
ex  ex  J  r 


(£)  -  fe) 


are  obtained  from  isentrcpic  flow  tables  as  a  function  of  jet  Mach  No.  and 
the  ratio  of  specific  heats  of  the  injected  gas. 

For  a  correctly  expanded  jet,  the  jet  exit  pressure  is  matched  to 
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the  air  pressure  immediately  ahead  of  the  jet,  and  the  jet  exit  diameter  is 
then  given  as 


d  =  5.  22 
ex 

where  T  -  540°R. 
o 

The  jet  exit  diameter  and  jet  Mach  No.  uniquely  determine  the 
maximum  penetration  distance.  Thus  the  relationship  obtained  from 
Figure  5  and  additional  calculations  is  shown  in  Figure  6  as 


(2.8) 


A 


=  dex  f<Mj’ 


M 


00 


=  90°) 


(2.9) 


The  change  in  penetration  distance  brought  about  by  a  reductier,  is  the 
injection  angle,  aQ,  is  shown  in  Figure  7.  Thus,  a  design  chart  can  be 
formed  (Figure  8)  using  Equation  2.  8  and  Figures  6  and  7  in  the  form  of 


=  a 


«Mj- 


M 


oo 


ao=90°) 


(2.  10) 


i 


2.2  Mass  Flow  Requirements 


Mass  flow  must  satisfy  the  following  three  requirements: 

1.  The  number  density  of  SF^  must  exceed  the  electron  number 

density  (xSF  »  X  el ). 

6 

2.  There  must  be  enough  SF^  to  obtain  a  fast  attachment  rate. 

3.  Enough  mass  must  be  injected  to  obtain  the  desired  penetration 
of  the  flow  field. 
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2,  3  System  Weight 


The  electrophilic  gas  is  stored  as  a  liquid  and,  in  most  applications, 
the  vapor  pressure  is  mor^  than  adequate  to  form  the  supersonic  gas  jet. 
Consequently,  the  injection  system  consists  of  a  pressure  vessel  for  gas 
storage,  nozzle,  on/off  valve,  and  metering  valves.  Since  small  pressure 
vessels  are  constructed  with  typical  weight-to-volume  ratios  of  0.05  lbs/in  , 
a  storage  vessel  to  hold  0.  5  lb.  of  liquid  SF^  (specific  gravity  =  1.6)  would 
weigh  1.  5  lbs.  The  valves  would  add  another  1.  5  lbs.  thus  giving  a  total 
system  weight  of  3.  0  lbs.  for  0.  5  lb.  of  SI^  or  a  system-to-gas  weight 
ratio  of  6.  0.  Specific  system  weights  will  be  given  in  Section  5  for  a 
Trailblazer  II  vehicle. 
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3.  NON- EQUILIBRIUM  COOLING 


In  order  to  find  a  means  of  removing  the  plasma  sheath  electrons 
blocking  communications  during  re-entry.,  the  thermodynamic  state  of 
the  ionized  plasma  layer  must  be  known,  and  then  a  way  must  be  found  to 
cool  this  flow  to  a  state  in  which  few  electrons  are  present.  The  flow  of 
air  over  the  body  surface  can  be  assumed  to  be  in  any  of  the  three 


following  states: 


1. 

Perfect  gas 

2. 

Equilibrium  gas 

3. 

Non-equilibrium  gas 

Although  the  perfect  gas  assumption  can  contribute  no  electrons  to  the 
flow,  it  is  nonetheless  useful  since  it  is  the  first  step  towards  more 
complex  assumptions. 

An  analysis  of  the  one-dimensional  mixing  and  cooling  of  a 
streamtube  with  the  assumption  of  both  a  perfect  and  an  equilibrium  gas 
has  been  completed  in  a  previous  report  (Ref.  1).  The  final  state  in 
the  cooling  process  for  the  perfect  gas  case  was  achieved  directly, 
although  the  final  temperature  had  to  be  specified  to  achieve  a  solution 
for  the  equilibrium  gas  case.  A  method  of  analysis  must  now  be 

f 

determined  for  the  more  realistic  condition  of  non-equilibrium  cooling. 
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3.1 


Two  Streamtabe  Formulation 


The  non- equilibrium  flow  of  a  Nitrogen-Oxygen  mixture  down  a 
streamtube  can  be  calculated  with  the  equations  presented  by  Vincenti 
(Ref.  5)  if  the  pressure  along  the  streamtube  is  specified.  The  formulation 
for  this  non- equilibrium  flow  has  been  programmed  on  a  7090  computer, 
and  thermodynamic  properties  as  a  function  of  streamtube  distance  can  be 
derived  by  following  the  various  chemical  reactions  in  the  streamtube. 
Cooling  the  flow  with  allowances  for  non- equilibrium  conditions  can  be 
determined  by  including  in  the  program  a  psuedo  reaction  describing  the 
interchange  of  energy  between  an  injected  coolant  and  non- equilibrium  flow. 
This  is  done  schematically  in  the  following  way: 


The  non- equilibrium  Nitrogen-Oxygen  reactants  are  contained  in  the  air 
streamtube  flowing  adjacent  to  the  coolant  streamtube.  The  coolant  and 
the  airstream  exchange  energy  in  the  mixing  region  until  they  both  reach  the 
same  final  energy  state;  however,  the  coolant  is  assumed  to  remain  a  non¬ 
reacting  gas.  Conservation  of  energy  and  momentum  between  the  two 
adjacent  streamtubes  is  stated  below:  ? 

<ina  +  CQdHc  =  0  (3.1) 

(3.2) 


du  +  C^du  =  0 
a  Q  c 
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Both  the  air  and  coolant  are  assumed  to  undergo  a  relaxation  process 
in  the  mixing  region  as  follows: 

“a  *  uf  +  (“a,  -  uf>  «"  ^ 
uc  =  uf(1.  0  -  e~Xx)  - 

where 

scaling  factor 
final  cooled  state 
streamtube  distance 

The  final  velocity,  uf,  is  unknown;  however,  by  forming  the  derivatives 
dua/ dx,  duc/dx 

and  substituting  these  values  in  Equation  3.2,  Uf  can  be  determined  to  be: 
u  \ 

ai  s 

TTT,  -  j  (3.6) 

1  / 


(3.4) 

(3.5) 


This  is  the  same  value  for  the  final  velocity  as  that  found  for  both  the  perfect 
and  equilibrium  gas  cases.  That  this  is  so  is  a  check  on  the  validity  of 
Reference  1  since,  regardless  of  what  non-isentropic  process  is  occurring 
between  two  stations,  the  momentum  equation  holds  considering  flow 
•  properties  at  stations  far  from  this  region. 
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The  relaxation  process  for  the  temperature  of  the  airstieam  and 
the  coolant  can  be  postulated  as  follows: 


Ta  =  Tf  +  (Tai_Tf)e"2XX  (3.7) 

Tc  =  Tf  '  (Tf-Tc  >  e’2XX  (3.8) 

ua=uf  +  (uai  -uf>e"XX  <3*4> 

uc  =  uf(l  -  e'Xx)  ■  (i.5) 


where  K  is  a  scaling  factor  which  determines  the  rate  of  cooling. 

Derivatives  of  du  /dx,  du  /dx  and  dh  /dx  are  now  formed  as  a  function  of 
a  c  a  * 

the  temperature  and  velocity  difference. 

Unlike  the  air  and  coolant  streamtube  velocities  the  temperature  has 

-2\x 

been  assumed  to  decay  as  e  since  this  assumption  yielded  the  simplest 
closed  form  solution  for  both  the  final  temperature  and  velocity.  The 
derivative  for  the  static  enthalpy  of  the  airstream  can  be  formed  using  the 
definition  for  the  static  enthalpy  of  a  calorically  perfect  gas, 


(3.11) 


and  eliminating  Tf  from  the  derivative  (see  Appendix  A).  The  result  is 
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In  the  same  manner,  derivatives  for  du  / dx,  and  du  / dx  can  be 

3.  C 

formed  from  Equations  3.4  and  3.  5  and  the  sum  of  these  two  equations. 
These  derivatives  are: 


du 

a 

"ET 


\C. 


Q, 


T+C 


Q, 


(u  -u  )' 
a  c' 


(3.  13) 


du  . 

~3T  =  +  1  +  Cn  (ua"uc) 


(3.  14) 


Having  formed  derivatives  of  static  enthalpy  for  air,  and  derivatives 

of  velocity  for  both  air  and  coolant  as  a  function  of  the  difference  in  coolant 

and  air  stream  temperatures  and  velocities,  a  numerical  iteration  procedure 

can  be  established  to  determine  new  values  of  air  and  coolant  velocities 

and  airstream  static  enthalpy.  The  static  enthalpies,  with  pressure  assumed 

constant  for  both  streamtubes,  will  yield  new  values  of  the  coolant  and  air 

temperatures  which,  in  turn,  will  allow  calculation  of  a  new  value  of  dh  / dx. 

a 

As  the  difference  (T  -T  )  and  (u  -u  )  approach  zero,  the  iteration  will  stop 

a  c  a  c 

and  these  values  will  be  the  final  ones. 

3>  2  Freezing  Criterion 

A  cooling  rate  matching  or  exceeding  the  chemical  reaction  rates  of 
the  air  species  comprising  lie  air  streamtube  must  be  determined  to  get  an 
equilibrium  cooling  solution.  Since  electrons  are  the  most  important 
chemical  specie  from  a  communication  standpoint,  a  freezing  criterion  will 
be  chosen  which  will  keep  electrons  in  equilibrium  while  cooling  is  occurring. 
To  do  this,  an  equilibrium  electron  concentration  and  its  rate  of  change  as 
a  function  of  distance  must  be  found  and  compared  to  the  rate  of  change  of 
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electrons  given  by  the  following  formula  taken  from  Appendix  B: 


=  nhh  [kf(N)(0)  -  Ve  )2J 


(3.15) 


where: 

(e  )  -  specie  in  particles/cc 

=  forward  reaction  constant  (cm^/mole  sec)  for 
N  +  O  —  NO+  +  e" 

*  backward  reaction  constant  (cm  /mole  sec)  for 
N  +  O  —NO+  +  e" 
u  =  air  streamtube  velocity. 

Since  the  primary  interest  is  in  cooling  cases  for  which 


^d  equilibrium  rate  of  change  of  electrons  will  be  compared  with 


(*). 


Ve‘)2 


(3.16) 


The  assumed  equilibrium  electron  concentration  is  found  in  the 
following  way.  The  final  cooled  equilibrium  state  is  known  from  previous 
research,  (Ref.l)  and  the  equilibrium  properties  of  the  streamtube  are 
expected  to  vary  in  a  way  similar  to  the  previously  obtained  perfect  gas 
cooling  solutions.  Therefore,  an  equilibrium  cooling  case,  in  which  the  final 
state  is  known  by  using  the  methods  presented  in  Reference  1,  can  be  compared 
with  a  perfect  gas  cooling  calculation  at  the  same  initial  enthalpy  (shown  in 
Figure  9).  Thus  the  assumed  variation  of  the  thermodynamic  variables  for 
the  equilibrium  case  can  be  determined.  The  enthalpy,  temperature,  pressure 
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and  velocity  thus  found  (Figure  10)  allows  calculating  electrons  as  a  function 
of  non- dimen sionalized  distance  \x  (Figure  11)  from  plasma  charts.  The 
equilibrium  electron  concentration  thus  found  is  then  differentiated  to 
determine 


where  \  *  represents,  in  effect,  the  distance  required  to  achieve  equilibrium 
recombination.  The  right  hand  side  of  Equation  3.16  can  now  be  calculated 
since  the  equilibrium  electron  variation  is  known  and  since,  k^,  the  backward 
rate  constant  in  the  reaction  describing  the  production  of  electrons,  is  simply 
a  function  of  temperature.  By  comparing  the  postulated  equilibrium  variation 
of  electrons 


(3.  18) 


with  the  actual  rate  of  change  of  electrons  as  defined  by 

de-  „  ye~)2 

dx  N  •  u 
o 

\  can  be  selected  such  that 


(3.19) 


I  (ie'\  <-  kb(<5")Z 

X.  y  dx  Je q  —  N  *u 


(3.20) 


and  thus  have  equilibrium  flow  cooling,  or  \  can  be  selected  such  that 


X  (xGt) 


eq 


and  thus  have  frozen 


Ve‘> 

>  "TTmT 


cooling. 


2 

(3.21) 

This  comparison  is  shown  in  Figure  12. 
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3.3 


Numerical  Results 


As  a  check  on  the  7090  computer  program,  runs  were  made  attempting 
to  duplicate  the  computations  made  on  a  Recomp  HI  computer  for  a  known 
perfect  gas  case.  This  was  done  by  setting  the  rate  constants  for  the  non¬ 
equilibrium  reactions  equal  to  zero  so  that  the  equations  describing  the  non¬ 
equilibrium  flow  reduce  to  a  perfect  gas  case  and  a  valid  comparison  could 
be  made.  The  results  of  these  runs  are  plotted  in  Figure  13  with  the  perfect 
gas  results  from  Reference  1.  The  initial  difference  in  the  temperature  and 
velocity  curves  is  due  to  a  difference  in  the  choice  of  computation  step  size 
for  the  7090  and  Recomp  III  programs. 

Similarly  a  check  was  made  for  an  equilibrium  gas  cooling  case  run 
at  the  condition  shown  in  Figure  14.  By  using  the  results  of  an  assumed 
equilibrium  variation  of  electrons,  derived  by  knowing  the  thermodynamic 
states  during  mixing  for  the  similar  (constant  initial  enthalpy)  perfect  gas 
case,  a  cooling  rate,  was  picked  such  that  the  flow  remained  in  equilibrium 
throughout.  With  this  \  =  .  01  cm  Figure  13  shows  that  the  final  cooled 
state  should  be  reached  at  approximately  \x  =  3.  0;  for  this  present  case  this 
would  mean  300  cm  downstream  of  initial  coolant  injection.  Since  the 

present  computer  program  involves  a  Runge-Kutta  integration  procedure 

-2 

which  requires  that  the  computation  step  size  be  <  1.10  cm  to  control 
the  propagation  of  errors,  only  the  first  1000°K  drop  in  temperature  was 
computed.  The  results  are  shown  in  Figure  14  and  compared  to  a  known 
equilibrium  solution. 

An  analysis  of  the  freezing  criterion  for  conditions 

T.  =  5075°K 

Pj  =  0.  1  atm 
CQ  =  .966046 

is  shown  in  Figure  15.  This  analysis  indicates  that,  unless  the  cooling 
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I 

i 

i 

I. 

.1  .  I 

length  \  is  greater  than  1000  cm,  non- equilibrium  cooling  would  occur. 

To  verify  this  observation,  computations  were  made  at  =10  with  d>e 

results  shown  in  Figures  16  and  17,  Although  the  temperature  dropped 

considerably  from  its  initial  value  (Figure  16)  the  various  specie 

concentration  (Figure  17)  remain  essentially  unchanged.  In  fact,  although 

2 

the  electron  concentration  dropped  by  a  factor  of  approximately  10  for 

the  non- equilibrium  case,  the  drop  would  have  been  much  greater  (i.e. , 

103}  had  equilibrium  occurred.  Therefore,  even  though  the  temperature 

is  cooled  within  10  cm,  the  electron  concentration  drops  only  by  a  factor 
2  5 

of  about  10  instead  of  the  10  reduction  possible  if  equilibrium  had  been 
reached.  A  cooling  length  of  1000  cm  would  be  required  to  achieve 
equilibrium.  This  is  an  unreasonable  distance  and  thus  shows  that  cooling 
techniques  are  not  an  effective  weapon  in  combatting  radio  blackout. 
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4.  ELECTRON  ATTACHMENT  TO  ELECTROPHILIC  GASES 


Essentially  electrophilic  gases  are  those  which  form  stable  negative 
ions  in  the  presence  of  electrons.  The  phenomenon  of  electron  attachment 
to  electrophilic  atoms  and  molecules  has  been  stud**  decades  and  must 

involve  energy  emission  since  the  negative  ion  is  at  ^  ,wer  energy  level  than 
the  original  neutral  species.  The  energy  emitted  to  attach  an  electron  is 
known  as  the  electron  affinity,  E^,  and  those  atoms  and  molecules  having  an 
electron  affinity  are  tabulated  in  Table  I;  chlorine  and  fluorine  atoms  "have 
the  highest  affinity  of  any  known  material. 

Electron  attachment  can  occur  either  by  three  body  attachment 
reaction  such  as 


F  +  X  +  e’  r  F'+X 

or  by  dissociative  attachment  reaction  such  as 


(4.1) 


SF6  +  e"  —  SF"  +  F 


(4.2) 


The  negative  ion  is  then  neutralized  in  collision  with  the  positive  ions. 
The  number  of  free  electrons  that  can  be  removed  from  a  plasma  can  be 
determined  from  equilibrium  calculations.  The  equilibrium  constant  is 
defined  as  the  ratio  of  the  concentration  of  products  over  the  concentration  of 
reactants.  For  three  body  attachment,  the  equilibrium  constant  is 


=  (F~) 
(F)  (e~) 


(4.3) 


and  for  two  body  dissociative  attachment,  the  constant  is 
(SF“  )  (F) 

K  =  - ? -  (4.  4) 

(SF6)(e-) 
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where  the  brackets  represent  molecules/cc.  The  equilibrium  constant  is 
related  to  the  partition  functions  by 

AE  " 

inKc  =  -TT  +  L  A7.inQ(A.)  (4.5) 

i=l 

where  A7^  is  the  difference  in  the  stoichiometric  coefficients  of  the  products 
and  the  reactants.  The  AEq  is  the  difference  in  the  zero  point  energy 
between  the  products  and  the  reactants.  The  partition  function  is  defined  as 
the  product  of  the  partition  functions  for  translation,  rotation,  vibration, 
and  electronic  energy  states,  and  is  given  below  for  translation  (Qfc)  and 
electronic  states  (Qe)>  respectively: 


(4.6) 


(4.  7) 


where  k  is  Boltzmann’s  constant* 

4.  1  Equilibrium  Attachment 

The  equilibrium  number  of  attached  electrons  is  independent  of  any 
attachment  reaction  and  is  a  function  only  of  temperature  and  pressure* 
Although  the  dissociative  attachment  reaction  is  a  faster  attachment  path 
than  molecular  dissociation  followed  by  atomic  attachment,  they  both  will 
eventually  reach  the  same  equilibrium  condition. 
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The  equilibrium  constant  for  the  three  body  fluorine  attachment 
reaction  can  be  written  using  Equations  4.  5  through  4.  7  and  the  constants 
in  Table  I,  as  follows: 


(4.8) 


The  difference  in  zero  point  energy  between  the  negative  ion  and  the  atom  is 
the  electron  affinity. 

The  electron  concentration  in  the  presence  of  halogen  atoms  was 
determined  from  the  simultaneous  solution  of  Equation  4.  8  and  the  following 
equations  for  conservation  of  atoms  and  charge: 


(F)  +  (F~)  -  (F)^  conservation  of  atoms 

(4.  9) 

(F~)  +  (e  )  =  (e  )  conservation  of  charge 


Note  that  (  )  denotes  initial  conditions. 

*  'o 

The  solution  of  these  equations  for  the  electron  concentration  is: 


The  usual  practice  is  to  have  the  fluorine  concentration  much  larger  than 
the  electron  concentration  since  under  this  circumstance,  fluorine  is  not 
appreciably  consumed  and  the  conservation  of  fluorine  atoms  can  be  neglected. 
Then  the  electron  concentration  is  given  as: 


(e“)  _  1 

r=r  ~  ic  m~  + 1 

(e  )  3cx  'o 

o 


(4.11) 


i 
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The  order  of  magnitude  of  electrons  removed  by  attachment  to  fluorine  is 
shown  in  Figure  18,  The  striking  feature  of  this  figure  is  that  almost  no 
electrons  are  removed  at  temperatures  above  3000°K.  The  results  are 
presented  in  terms  of  the  partial  pressure  of  fluorine;  the  mole  fraction 
of  fluorine  added  is  given  by: 

X  =  P'(F)/Po 

The  inverse  equilibrium  constant  for  other  electrophilic  gases  is 
shown  in  Figure  19.  The  halogen  atoms  are  shown  to  be  far  superior  to 
other  atoms  in  their  ability  to  attach  electrons. 


4. 2  Relative  Importance  of  Molecular  and  Electron  Collisions 
in  tHssociation  of  SF^  ' 


The  mode  of  electron  attachment  is  determined  by  the  relative 
speeds  of  molecular  and  electron  collisions  in  dissociating  SF^  which  occurs 
by  the  following  reaction: 


sf6  +  X  ~  sf5  +  F  +  X 

The  dissociation  rate  is  specified  for  Equation  4. 12  as 
d(SF,) 

TT  =  "  kD  <SF6><x> 


(4.  12) 


(4.  13) 


Integration  of  Equation  4. 13,  assuming  constant  air  temperature  and  density, 
gives  an  exponential  change  in  SF^  concentration,  e.g. : 


<SF6)  =  <SF6)o  e 


■kD<X>‘o 


(4. 14) 


where  (X)  is  the  air  concentration.  The  dissociation  rate  coefficient,  kjy 
is  obtained  from  experiments.  In  the  shock  tube  experiments  (Ref.  6), 


20 


\;C  63-78-R2 


complete  dissociation  of  SF^  was  observed  to  occur  at  some  time  prior  to 
the  ionization  of  air.  If  the  rate  determining  step  for  the  complete 
dissociation  of  SF^  is  the  first  step  (Equation  4. 12),  then  the  rate 
coefficient  can  be  specified  from  Equation  4,  14,  as  follows: 


ko  =7X5*7 


<SF6>o 

Tsftt 


(4.15) 


Substituting  the  observed  ionization  for  the  dissociation  time,  tj y  and 
assuming  that  the  SF^  is  99.  9  percent  dissociated,  the  dissociation 
coefficient  at  test  conditions  of  3840°K  and  1. 5  atm  is  given  as 


k^  =  4.3x10 


3  ,  ,  -1  -1 
cm  molecules  s 


The  temperature  dependence  of  the  rate  coefficient  can  be 
determined  by  assuming  a  simple  Arrhenius  law,  as  follows: 


k  -Q  —  A  e 


■  E/RT 


(4. 16) 


The  activation  energy  for  dissociation  is  the  energy  of  the  bond  broken  in 
the  reaction.  The  first  bond  is  the  strongest,  thus  substantiating  the 
assumption  that  the  first  step  determines  rate.  With  the  bond  energy  for 
Equation  4.  12  given  in  Table  II  as  85  kcal/mole,  the  rate  coefficient  can 
be  specified  as 


kp  =  3.  5  x  10‘*°  e-8500°/RT 


(4.17) 


which  is  within  a  factor  of  two  of  the  value  obtained  by  direct  observation  of 
Cl^  dissociation  (Ref.  7).  Oxidation  reactions  with  SF^  are  possible  and 
would  speed  up  the  dissociation  (Ref.  8). 

The  rate  coefficients  for  the  dissociative  attachment  reaction, 
Equation  4.2,  can  be  specified  from  electron  beam  experiments  used  to 
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measure  attachment  cross  sections.  The  cross  section  for  SF^  is  given  in 
Reference  9  and  is  a  function  of  the  electron  velocity.  The  rate  coefficient 
is  defined  as 


*2  = 


U  a(U)  dU  . 


(4.18) 


There  is  no  exponential  dependence  since  the  activation  energy  for  the 
dissociative  attachment  is  essentially  zero.  This  integral  can  be  evaluated 
as  a  function  of  temperature  assuming  an  appropriate  velocity  distribution. 
For  illustration  purposes,  approximating  the  integral  will  be  sufficient 
and,  thus,  the  rate  coefficient  can  be  derived  as  follows: 


H  « 


peak 


(4. 19) 


where 


tT: — n*L 

^  ’’’"’el 


el 

\ 

apeak  (^6^  =  2  x  10”^  cm^  @  E  =:  0. 1  eV.  (see  note) 

\ 

\ 

E  =  electron  energy  at  cross  section  peak. 

Thus  the  dissociative  attachment  rate  coefficient  is  given  as 


-  •  o  ]  Y 

k^  =  1.3xl0"li  cm  molecules  ~  s 


(4.  20) 


Loss  of  electrons  by  resonance  electron  capture  is  insignificant  at 
temperatures  above  1000°K  since  the  resonance  capture  cross  section 
consists  of  a  very  narrow  peak  centered  at  about  zero  electron  energy. 


Note  added  in  proof: 

Cross  sections  as  high  as  2.6  x  10"^  cm^  have  been  reported  by  Mahan 
and  Young  (Ref.  16).  This  would  greatly  increase  the  rate  of  dissociative 
attachment. 
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The  relative  importance  of  the  two  dissociative  processes  is  obtained 
by  plotting  the  ratio  of  the  molecular  dissociation  rate  to  the  dissociative 
attachment  rate: 


d(SF6) 

*  _  ~  & 

molecular 

3TSF-) 

clt 

electron 

-ysr6)(x) 

-ki,(SF6)(e”) 


(4.21) 


Substituting  for  Equations  4.13  and  4. 19, 


t  -  35  T'1/2  e-«000/T  X-1 

e 

the  mole  fraction  of  electrons  is  defined 


(4.22) 


X 

e 


_  (e“) 

"  T xT* 


(4.23) 


Equation  (4,22)  is  plotted  as  a  function  of  temperature  in  Figure  20 
to  illustrate  that  dissociative  attachment  is  the  dominant  reaction  for  typical 
values  of  the  electron  mole  fraction  at  temperatures  below  about  4000°K. 
Thus,  in  expanded  flow  regions  where  the  temperature  is  low,  the  SF^  will 
remove  electrons  by  dissociative  attachment.  In  the  high  temperature 
regions  near  the  stagnation  point,  molecular  dissociation  is  the  predominant 
reaction  and  electrons  can  or'  y  be  removed  by  three  body  attachment. 


4.  3  Dissociative  Attachment  Rate  in  Expanded  Flow  Regions 

On  many  vehicles,  the  communications  antennas  are  located  aft 
of  the  vehicle  nose  in  an  expanded  flow  region.  The  air  around  the  antenna 
will  have  expanded  from  a  stagnation  point  condition  and  will  be  cool; 
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howsver  the  electron  density  will  be  large  because  the  recombination  rates 
are  slow  compared  to  the  expansion  rate.  Thus,  a  large  electron  mole 
fraction  exists  in  a  cooled  gas.  An  electrophilic  gas  injected  into  th.'s  type 
of  flow  will  remove  electrons  by  the  dissociative  process  described  earlier. 

The  rate  of  electron  removal  by  dissociative  attachment  is  given  as 


d(e”) 

“3t 


=  -  4  (SF6)  (e“) 


(4.  24) 


_3 

where  (  )  represents  molecules  cm  .  This  expression  is  valid  only  until 
equilibrium  is  approached.  At  near  equilibrium,  reverse  reactions  will 
begin  to  be  significant.  The  amount  of  electron  density  reduction  is  obtained 
by  integrating  Equation  4.24  assuming  that  the  air  temperature  and  density 
remain  fixed.  The  concentration  of  SF^  is  also  unchanged  since  the  amount 
added  greatly  exceeds  the  electron  density.  Integrating  Equation  4.24  also 
yields  an  exponential  electron  density  decay  of  the  form 


(e")  =  (e")Q  exp  £ -k^(SF6)(At)  1 


(4.25) 


By  defining  the  concentration  of  SF^  in  terms  of  mole  fraction  of  the  air 
density  and  substituting  the  rate  coefficient  in  Equation  4.20,  the  elapsed 

time  to  reduce  the  electron  density  by  100  by  dissociative  attachment, 

*100  . 

,  is  given  as 


joo  v  ;  i,7xio“ 

2  3F6  = 


(4.26) 


The  attachment  time,  given  in  Equation  4.26,  is  plotted  in  Figure  21 

as  a  function  of  temperature  and  pressure  (the  Mollier  diagram  in 
Reference  10  was  used  to  obtain  density). 
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Three  Body  Electron  Attachment  Rate 


The  attachment  time  for  the  three  body  reaction  will  be  the  sum  of 
the  dissociation  time  of  the  halogen  compound  and  the  direct  attachment 
time  of  Equations  4.  12  and  4. 1.  The  elapsed  time  to  reduce  the  electron 
concentration  to  l/lOO  its  original  value  is  obtained  from  Equation  4. 14 
with  the  second  body  concentration  expressed  as  density.  The  expression 
for  t^^  is  given  below  and  the  values  are  plotted  in  Figure  22: 


.100  ,  ,n-12  -1  8600/RT 

tp  =2x10  p  e  ' 


(4.27) 


The  rate  of  three  body  electron  attachment  to  fluorine  is  defined  from 
Equation  4.  1  as 


~  ~  kj(F)  (e  )  (X) 


(4.28) 


The  rate  coefficient,  obtained  from  shock  time  measurements  by  Good 
(Ref.  6),  is 

k3  =  1.5  x  10‘37  e^4,  000/RT  cc2_molecule-2_g- *  (4.29) 

Integrating  Equation  4,28  yields  the  exponential  decay  of  electron 
concentration,  as  follows: 


(e~)  =  (e")Q  exp  f-k^F)  (X)  t3J 


(4.30) 
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The  time  required  to  reduce  the  electron  concentration  to  l/lOO  its  original 
value,  defined  as  ,  is  given  as 


100  _  (e") 

3  =  F3(F)(Xy 


7  x  10"®  e-«000A 


(4.31) 


The  attachment  time,  t^  f  given  in  Equation  4.  31  has  been  plotted  in 
Figure  23  as  functions  of  temperature  and  pressure. 

The  total  attachment  time,  representing  the  sum  of  the  dissociation 
and  three  body  attachment  time,  is  given  in  Figure  24.  At  conditions  where 
the  dissociation  time  is  short,  the  attachment  time  is  long  and  conversely. 
Consequently,  the  combined  time  is  always  long,  and  the  minimum  time  is 
reached  at  a  temperature  of  3000  K. 

The  total  attachment  time  is  never  less  than  10  ms  for  typical 
re-entry  pressures,  and  consequently,  three  body  atomic-electron 
attachment  is  a  mechanism  that  can  only  be  used  in  stagnation  point 
regimes.  However,  its  use  even  at  this  location  is  doubtful  since  the  weight 
requirements  to  cool  the  flow  to  3000°K  may  be  prohibitive. 
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5.  INJECTION  SYSTEM  APPLIED  TO  TRAILBLAZER  VEHICLE 


The  injection  system  concept  will  now  be  carried  to  specifics  and 
applied  to  the  Trailblazer  II  vehicle..  For  purposes  of  this  example,  the 
vehicle's  trajectory  will  be  considered  a  constant  velocity  flight  of 
18,  000  fps  between  125,  000  ft.  and  300,  000  ft.  Essentially  the  vehicle  is 
a  blunted  sphere-cone  having  a  6-in.  nose  radius,  a  9°  cone  angle,  and 
a  total  vehicle  length  of  26  in.  (Figure  25).  The  antenna  is  located  near  the 
rear  of  the  vehicle. 

The  plasma  on  the  Trailblazer  II  consists  of  a  thin  layer  with  a 
thickness  of  about  one  tenth  the  nose  diameter  (Ref.  11).  This  represents 
the  stagnation  point  flow  that  expands  around  the  vehicle  nose.  The  flow 
field  properties,  shown  in  Table  HI  as  a  function  of  altitude,  were  calculated 
by  assuming  equilibrium  expansion  of  the  flow  up  to  the  sonic  line. 
Thereafter,  the  flow  was  assumed  to  be  frozen  at  the  sonic  conditions  and 
expanded  according  tc  expressions  given  by  Heims  (Ref.  12).  The  pressure 
is  assumed  to  expand  to  one  tenth  the  normal  shock  pressure  as  indicated 
by  Yakura  (Ref.  13)  for  hypersonic  flow  over  blunted  cones. 

5.  1  Mass  Flow  Injection  Rates 

The  attachment  reaction  order  and  speed  will  depend  on  whether  the 
electrophilic  gas  is  injected  into  a  high  or  low  temperature  flow.  If  the 
electrophilic  gas  is  injected  into  the  stagnation  point  region,  the  high 
temperatures  will  rapidly  dissociate  any  existing  electrophilic  molecules 
thus  leaving  the  atoms  to  attach  electrons  by  a  three  body  reaction.  If 
the  electrophilic  gas  injected  into  an  expanded  flow  region  where  the  high 
temperatures  have  cooled,  the  electrophilic  molecules  will  remain  intact 
and  will  be  free  to  reduce  the  plasma  frequency  by  dissociative  electron 
attachment. 
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5.  1. 1  Three  Body  Attachment 


The  three  body  attachment  process  is  normally  quite  slow  and  not 
useful  in  alleviating  the  plasma  blackout;  however,  at  high  pressures  the 
process  is  useful.  In  general,  the  allowable  attachment  time  is  fixed  by 
the  flow  time  of  the  electrophilic  gas  between  the  point  of  injection  and  the 
antenna.  Thus,  the  attachment  time  is  defined  as 


t  < 


j.2.5 


ds 

IT 


(5.1) 


where  ds  is  the  distance  increment  along  the  body.  The  Trailblazer  II  body 
distance  from  stagnation  point  to  the  rear  antennas  is  2.  5  ft. 

For  three  body  attachment,  the  attachment  time  required  to  reduce 
the  electron  density  to  one  percent  of  its  original  value  is  defined  as 


100 

h 


(5.2) 


where  (t^^  X^)  is  the  modified  attachment  time  shown  in  Figure  23. 
Substituting  Equation  5.2  into  5.  1  results  in  the  following  attachment  time 
criterion: 


(t 


100 

3 


V  <  yF 


ds 

TT 


(5.3) 


The  integral  has  the  value  of  2.  8  x  10  seconds,  so  that 


(t^00  Xp)  <  2.8  x  10“4  Xp 

Considering  for  the  moment  that  the  upper  limit  of  material  addition  is 
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=  0.  1,  the  modined  three  body  attachment  time  is 
(t*°°  Xp)  <  2.8xl0'5 

and  as  indicated  in  Figure  23,  this  is  achieved  only  at  pressures  above 
10  atm  even  with  the  low  temperature  of  T  =  2000°K.  For  the 
TraiJblazer  trajectory,  three  body  attachment  will  only  be  possible  at 
altitudes  below  175,  000  ft. 

5.1.2  Two  Body  Attachment 


The  only  reaction  which  will  remove  electrons  during  the  entire 
re-entry  trajectory  is  two  body  dissociative  attachment,  and  SF^  is  the 
most  effective  gas  for  this  since  its  attachment  cross  section  is  the  largest 
known  and  its  equilibrium  number  of  electrons  that  can  be  attached  is  also 
large.  The  dissociative  attachment  will  reduce  the  electron  density  to  at 
least  the  equilibrium  concentration  level  of  pure  air. 

To  avoid  dissociation  of  SF^  when  it  is  injected  into  the  plasma, 
the  flow  field  temperature  at  the  injection  point  must  be  less  than  3000°K. 
For  all  altitudes  above  150,  000  feet,  this  can  be  achieved  by  injecting  at 
the  shoulder;  however,  at  lower  altitudes,  some  SF^  dissociation  must  be 
permitted  to  cool  the  plasma  so  the  remaining  SF^  will  attach  electrons. 
Details  such  as  this  remain  to  be  solved  during  an  actual  design.  The 
calculations  presented  below  only  illustrate  system  feasibility. 

The  amount  of  SF^  to  be  added  is  determined  by  three  criteria,  as 
mentioned  earlier:  electron  density,  attachment  rate,  and  jet  momentum. 
The  SF^  concentration  added  should  be  greater  than  the  electron 
concentration,  and  for  design  purposes,  the  mole  fraction  of  SF^  added  to 
the  plasma  should  equal  100  times  the  electron  mole  fraction,  i.  e. 


=  100X 

e 


(5.4) 
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The  mass  flow  rate  is  specified  as 


ITl  gjp 

Ai  =  100  X  ,  p  U  Ad  - — 

el  rs  s  P  m  - 
air 


(5.5) 


where 


U 


m 


SFZ 


m 


air 


3 

is  the  density/,  Ibs/ft 
is  the  sonic  flow  velocity 
is  the  molecular  weight  of  SF^ 
is  the  molecular  weight  of  air 

is  the  streamtube  area  containing  the  plasma  passing 
over  the  antenna 


The  effective  plasma  thickness  on  a  9°  blunted  cone  is  taken  as  the 
maximum  value  tabulated  by  Russo  (Ref.  11)  in  his  Table  I  for  a  9°  cone 
at  M  =  20  between  60,  000  and  230,  000  ft.  altitude. 


dp/Dn  <  °«i44  (5.6) 

where  the  nose  diameter,  Dr,  for  the  Trailblazer  His  given  as  12  inches. 
Thus  the  plasma  area,  Ap,  is  defined  as  the  streamtube  cross  section 
containing  the  plasma  which  affects  the  antenna.  The  area  is  the  plasma 
thickness  dp  and  the  required  transmission  window  width  for  the 
Trailblazer  II.  The  cross  sectional  area  at  the  antenna  will  thus  be 
essentially 

2 

Ap  =  dp  x  w  «  6  inches  ^ 

The  ?  ea  of  the  same  streamtube  at  the  point  of  injection  at  the  shoulder  will 
be  three  square  inches.  The  -equired  mass  flow  rate  is  shown,  in  Figure  26. 
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Thus  the  electron  density  mass  flow  rate  requirement  is  given  as 


m  =  1 . 2  ::  10“ 


PSXel  =  ‘-24 


x  10"20  N 


el 


(5.8) 


where  Ns  is  the  sonic  line  number  density. 

The  amount  of  SF^  required  for  dissociative  attachment,  of  99 
percent  of  the  original  electrons  within  a  2  ft.  length  is  specified  as 


100 

*2 


< 


ds 

"V 


(5.9) 


and  substituting  for  the  time  with  Equation  5.  2 


<  2.  8  x  10 


-4 


and 


3.  5  x  10' 


SF6 


) 


(  5.10) 


(5.11) 


The  mass  flow  rate  of  SF^  that  must  be  injected  according  to  the  attachment 
rate  requirement  specified  by  Equation  5.  11  is  given  as 


m 


SF, 


m  =  p  U  A 
Ks 


P  m 


air 


[3-5xl°3  (  ^°°  XsfJ] 


(5.12) 


and  substituting  the  previous  values  of  Ap  andl)  and  defining  weight  density, 
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pg,  in  terms  of  the  number  density,  Ng, 


m  =  2.2  x  10"26  Ng  ^2°°  Xsf6)  lbs/s  (5. 13) 

For  the  Trailblazer  il  conditions  given  in  Table  III  and  Figure  21,  the  product 

N  (tl^x  )  never  exceeds  1  x  10^  molecules- s  ^  Thus  the  flow  rate 
s  2  SF /  r,  3 

®  ^  .o 

required  for  attachment  is  always  less  than  10  7  lbs/s  which  is  negligible 

compared  to  the  electron  density  requirement.  Furthermore,  the  reduction 

of  the  electron  density  to  0. 1  percent  of  its  original  value  would  require 

only  a  50  percent  increase  in  time  which  is  still  negligible  providing 

equilibrium  is  not  approached. 

On  the  basis  of  momentum  criteria,  the  mass  flow  requirements  are 
obtained  from  Figure  8.  The  jet  centerline  penetration  distance  is  defined 
to  be  half  the  plasma  layer  thickness. 

A  =  dp  =  0.9  inches 

The  mass  flow  required  to  penetrate  0.9  in.  into  the  flow  at  the  shoulder 
(M  =:  2.  0)  using  a  jet  Mach  No. ,  M.  =  2,  is  found  to  be 

m.  =  1.  1  x  10-2  Pg  lbs/s  {5.14) 

where  the  jet  exit  pressure  is  matched  to  the  local  shoulder  pressure  P  , 

5.2  Reduction  of  Plasma  Attenuation  and  Reflection 


Since  mass  flow  rates  are  fixed  by  the  jet  momentum  requirements, 
there  will  be  more  than  sufficient  time  and  SF^  concentration  to  bring 
the  plasma  sheath  to  equilibrium  air  conditions  within  a  distance  of  two  feet. 
Figure  27  illustrates  the  dramatic  alleviation  of  S-band  blackout  by  the 
addition  of  SF^.  The  signal  attenuation  is  reduced  to  zero  and  the  reflection 
losses  reduced  to  a  3  db  peak  at  130,  000  ft. 
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5«  3  Gas  Supply  System 


The  required  mass  flow  rates  are  shown  in  Figure  26.  The 
momentum  requirements  of  the  Mach  2  jet  predominate  and  determine 
the  mass  flow  for  the  Trailblazer  II.  The  total  mass  required  is 
determined  as 


S  125, 000 

m.  =  0.28  lbs  (5.  15) 

300,000  3  °° 

The  SF^  can  be  conveniently  stored  under  its  own  vapor  pressure  of 
310  psig.  and  the  amount  of  gas  that  must  be  stored  is  determined  by  the 
requirement  that  following  expulsion  of  0.  28  lbs  of  gas  at  125,  000  ft,  the 
storage  vessel  must  still  have  a  pressure  of  Pgx  =  Pg  =  84  psi.  The  amount 
of  gas  required  for  isentropic  expansion  is  determined  from  the  following 
relations: 

p  =  (5. 16) 


where 

m  is  the  weight  of  gas  in  storage 
V  is  the  storage  volume 
T  is  the  gas  temperature 
R  is  the  gas  constant 
P  is  the  gas  pressure 
Differentiating  Equation  5.16, 

dm  _  V 

7KT  =  IT 


(5.17) 
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The  differential  is  evaluated  from  the  isentropic  relation 


lli 


«)■(*) 


(5.  18) 


and  substituted  into  Equation  5. 17  to  obtain 


dm  _  m 

"3t“  =  ~P 


2(7-0 


7 


dP 

Si 


which  is  integrated  as 


(5.19) 


(5.20) 


.With  the  requirements  that  Pq  equal  300  psia,  Pf  equal  84  psia,  and 
7  equal  1. 09,  the  ratio  of  final  to  initial  mass  weight  of  SF^  stored  in  the 
vessel  is  found  to  be 


and  the  expelled  weight  is  defined  as 
-  m^  -  0.  28  lbs  . 

Thus  the  initial  stored  weight  is  determined  to  be 
m  =  0.  4  lbs. 

G 
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The  0.4  lbs  of  SF^  required  for  the  Trailblaser  II  flight  test  is  a 
representative  value.  Optimization  is  required  on  the  orifice  size,  jet 
Mach  No.  and  gas  weights  to  cover  the  wide  range  of  exit  pressures 
required  for  a  correctly  expanded  jet. 

The  density  of  the  SF^  stored  as  a  liquid  under  its  own  vapor 

pressure  at  70°F,  is  3.25  lbs/ft^,  and  the  storage  vessel  size  required 

3  3 

for  0.4  lb  of  SF^  is  0.  124  ft  or  215  in  with  a  3.7  in.  radius.  A  spherical 
vessel  this  size  made  of  4340  steel  with  a  burst  pressure  of  900  psi 
would  weigh  1.5  lbs,  and  the  regulator,  fill  valves,  on-off  switches, 
squib  valves,  and  required  plumbing  would  add  another  1.5  lbs  to  the  system 
(Ref.  15).  Thus  the  total  system  would  weigh  3.4  lbs  when  charged  with 


The  system  parameters  that  have  been  specified  illustrate  the 
feasibility  of  alleviating  radio  communication  blackout.  This  is  by  no 
means  the  optimum  or  desired  system  for  a  flight  test.  The  specific 
program  may  have  smaller  volume  available.  In  which  case,  the  amount 
of  mass  and  storage  volume  can  be  reduced  by  pulsing  the  gas  jet. 

5.4  Jet  Effect  on  Vehicle  Aerodynamics 

The  use  of  a  side  jet  to  create  an  electromagnetic  window  will  create  a 
moment,  upsetting  the  vehicle  trim.  The  trim  angle  can  be  obtained  by 

balancing  the  jet  moment  with  the  aerodynamic  moment  produced  by  the  trim 
angle. 


/ 

/ 


V 


=  c 


..2 
P  V 
00  00 


7 rD 


N 


a 


0 '  .  X 

trim 
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where 


N.  is  the  jet  normal  force 
J 

i  is  the  axial  distance  between  the  jet 

and  the  center-of-pressure. 

D  is  the  vehicle  base  diameter 


is  the  normal 


or  .  is  the  trim  angle 
trim  & 


is  the  static  margin 


For  the  TrailbAazer  II  vehicle 


=  5.  4  inches 


=  .015  per  degree 


20  inches 


=  0.  54  inches 


The  jet  normal  force  is  assumed  to  be  twice  the  jet  exit  momentum 


N.  =  2  m  U. 
J  J 


For  the  Mach  2  jet  assuming  a  gas  storage  temperature  of  520  R, 

N.  =  50.5  m  lbs 
J 

Thus  the  trim  angle  can  be  evaluated  as  a  function  of  altitude  using  the 
following  relation 


a*  .  =214 

trim 


P  V 
00  00 
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For  the  mass  flow  rates  specified  in  Figure  27,  the  trim  angle  is  calcu¬ 
lated  to  be  about  0.35  degrees  down  to  180  km  and  hereafter  increases  to 
0.6  degrees  at  125  km. 

Such  a  value  for  the  trim  angle  of  attack  is  not  considered  a 
problem. 


6.  CONCLUSIONS 


The  use  of  a  high  heat  capacity  gas  to  cool  the  plasma  and  create  a 
radio  communications  window  was  shown  to  have  no  effect  on  the  plasma 
frequency  or  electron  density.  The  use  of  such  a  gas,  however,  does  cool 
the  plasma  temperature  but  has  no  effect  on  the  electron  recombination 
rates  which  remain  too  slow  to  appreciably  reduce  the  electron  density 
near  the  point  of  gas  injection.  Numerical  calculations  at  T  =  5075°K  and 
P  =  0.  1  atm  indicate  that  equilibrium  cooling  of  electrons  would  require  a 
distance  of  about  ten  meters. 

Speedy  removal  of  electrons  can  only  be  achieved  by  adding 
material  to  react  with  the  electrons  to  form  new,  fast  recombination 
mechanisms.  Halogens  are  particularly  useful  for  this  and  the  electrophilic 
gas,  SF^,  has  the  largest  known  electron  attachment  cross  section.  Thus, 
SF^  added  to  the  plasma  will  rapidly  remove  electrons  in  a  two  body 
process  known  as  dissociative  attachment.  Studies  to  determine  conditions 
where  dissociative  attachment  would  occur  showed  that  if  the  SF^  is 
injected  into  a  plasma  which  has  cooled  by  expansion  to  temperatures  below 
3000°K,  the  SF^  molecule  will  remain  intact  until  it  encounters  a  free 
electron.  Typical  attachment  times  for  re-entry  plasma  conditions  are 
less  than  ten  microseconds  for  ten  percent  mole  fraction  of  SF^  added. 

A  SF^  gas  injection  system  for  the  Trailblazer  II  has  been  studied 
and  demonstrated  feasible.  The  preliminary  design  developed  during  this 
study  indicates  that  a  system  weighing  3.4  pounds  is  all  that  is  required  to 
reduce  the  S-band  attenuation  and  reflection  to  essentially  zero  throughout 
the  re-entry.  A  detailed  design  is  now  required  to  develop  the  actual 
mechanical  design  of  the  storage  vessel,  regulator,  nozzle  and  associated 
plumbing.  For  the  Trailblazer  II  vehicle  it  appears  feasible  that  the  SF^ 
can  be  stored  as  a  liquid  under  its  own  vapor  pressure  in  tubing  wrapped 
around  inside  the  vehicle  base.  Thus  the  gas  injection  system  can  be 
designed  and  installed  which  will  not  affect  the  payload  carried. 
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TABLE  I 


PROPERTIES  OF  ELECTROPHILIC  GASES 


Species 

Energy 

Level 

n 

Electronic 

Degeneracy 

gn 

Electronic 

Energy 

V  ev 

Ionization 

Energy 

ev 

Electron 
Affinity 
ev  i 

F 

0 

4 

0 

17.418 

-3.  63 

1 

2 

12.71 

Cl 

0 

4 

0 

13.01 

-3.78 

1 

2 

0.109 

2 

1 

8.825 

Br 

0 

4 

0 

1 1 . 84 

-3.  54 

I 

0 

4 

0 

10.454 

-3.  24 

O 

1 

5 

0 

13. 614 

-1.0 

2 

3 

1.967 

SF6 

0 

CF3 

10.  1 

-  .9 

h2° 

9.511 

s 

10.357 

-2.  5 

c 

11. 256 

-1.  13 

41 


MC  63-78-R2 


TABLE  II 


PROPERTIES  OF  SF6  AND  CF4 


Specie 

Molecular 
Weight 
g  /  mole 

M 

Specific 

Heat 

cal/mole-  c 
CP 

Enthalpy 
at  3006K 
cal/mole 

H300 

Bond 

Energy 

k/cal/mole 

AE 

AH 

Dissociation 

kcal/mole 

SF6 

146.06 

23.  22 

6945 

85 

,4  26 

/ 

SFS 

127.06 

4R 

— 

65 

— 

SF4 

108.06 

4R 

— 

65 

— 

SF3 

89.06 

4R 

— 

65 

SF2 

70.06 

4R 

— 

65 

— 

SF 

51.06 

7/2R 

— 

65 

— 

S 

3  2.06 

5.66 

— 

-- 

j  F 

.  19 

5.436 

— 

1 

— 

CF4 

88.01 

4R 

5992 

123 

462 

CF3 

69.01 

4R  ! 

113 

— 

CF2 

50.01 

4R 

— 

113 

— 

CF 

31.01 

7/2 

— 

113 

— 

C 

\  1  2.  0 1 

4.98 

— 

— 
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Figure 


x/Dn 


1.  Radio  Attenuation  Around  a  Blunt-Body  Vehicle 
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REGION  OF 
MIXED  FLOW 


Figure  3.  Plan  View  of  Jet  Injection  System 
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TRAILING  EDGE 
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Figure  5.  Centerline  Penetration  Distances 


now  (A/«) 


Figure  6.  Effect  of  Jet  Mach  Number  on  Penetration  Distance 


A/Aa  =90° 


Figure  7 


Effect  of  Jet  Injection  Angle  on  Penetration  Distance 


\ 
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Figure  9.  Mass  Flow  Requirements  to  Cool  Plasma  to  2400  K 


<0  (0  « 

6  6  <5 

SlINfl  AUVUJ.I8UV 
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Figure  10.  Estimate  of  Equilibrium  Gas  Variation  of  Temperature,  Pressure,  and  Enthalpy 


(ELECTRONS  CM* 


Figure  13.  Perfect  Gas  Cooling 


5: 


Equilibrium  Gas  Cooling 


60 


MC  63-78-R2 


igure  17.  Specie  Concentration  During  Frozen  Flow  Cooli 


ELECTRON  DENSITY  WITH  FLUORINE 
ELECTRON  DENSITY  WITHOUT  FLUORINE 


MOLECULAR  DISSOCIATION  RATE 


Figure  20.  Relative  Importance  of  Molecular  and  Electron 
Collisions  in  SF^  Dissociation 
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TEMPERATURE* 4000  °K 
3000  °K 
2000  °K 


PRESSURE  (ATM) 


Figure  21.  Dissociative  Electron  Attachment  Time 
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Xr  FLUORINE  MOLE  FRACTION 


'ATTACHMENT  TIME  TO  REDUCE 
ELECTRON  DENSITY  BY  IOO 


PRESSURE  (ATM) 


TEMPERATURE 


Figure  23.  Electron  Attachment  Time 
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S0NO03S  4  X  + 


o  a 


PRESSURE  (ATM) 


Figure  24.  Total  Time  to  Form  Atomic  Fluorine  and  Attach  Electrons 
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Figure  27.  S-Band  Transmission  Losses 


APPENDIX  A 


The  derivatives  dh  /dx,  du  /ax,  du  /dx  to  be  derived  will  be  exact 

cl  cL  C 

for  a  perfect  gas  since  they  were  derived  with  this  assumption;  however, 
the  form  of  the  derivatives  will  be  changed  so  that  they  can  be  used  regard¬ 
less  of  the  thermodynamic  state  of  the  hot  gas  to  be  cooled.  To  form 
dha/dx,  the  derivative  is 


dh 

a 

Hx" 


dT 

.  a 
dx 


(Al) 


Differentiating  Equation  3.7  gives 


(A2) 


This  derivative  must  now  be  put  in  a  form  which  is  a  function  of  the 
temperature  difference  of  the  air  and  coolant  streamtube.  In  this  way,  the 
derivative  will  have  been  defined  with  the  assumption  of  a  perfect  gas  and 
will  still  be  applicable  to  any  thermodynamic  state  as  long  as  the  function 
(T  -T  )  is  being  continuously  determined.  To  eliminate  (T  -T.)  from 

^  C  2>  X 

Equation  A2,  the  perfect  gas  energy  evaluatior  of  Equation  3. 1  is  used. 
First,  this  equation  is  rearranged  to  get, 


+  C 


u 


1 
T 


mr. 


(A3) 
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Multiplying  this  equation  by  -1.0  and  adding  (C0  +  Cn  C—  )  T  to  both 

a  Ul  c  al 

sides,  it  can  be  solved  for  (T  -T-),  and  the  following  aquation  can  be 

2i .  1 

obtained:  1 


2  C, 


(vT'h 


QiCpc(T»rTc»)'  ^  TTikl 


^Pa  -'Q,  VP„ 

3.  I  C 


Substitution  of  Equation  A4  into  Equation  A2  gives: 


2\  Cp  e 
a 


-2\x 


2 

CQiCpc(T*rTci)'^ 


1  +  c 


a 


“  2  \x 

which  is  still  not  in  the  form  desired,  but  by  substituting  for  e  its 

equivalent  as  a  function  of  (T  -T  )  the  final  form  of  dh  /dx  will  be 

a  c  «  a' 

determined.  If  Equation  3.7  is  subtracted  from  3.8 


(Ta-Tc)  = 


/t  -T  \ 
{  al  Cl) 


-2\x 


-2\x 


Substituting  the  value  of  e  from  the  above  equation  into  Equation  A5 

gives,  after  rearrangement,  the  final  form  of  dh  /dx. 

a  ’ 


2X.  C, 


1  +  C, 


al  C1 


(Ta”Tc) 


!1 


72 


MC  63-78-R2 


The  derivatives  of  air  and  coolant  streamtube  velocity  with  stream- 
tube  distance  can  be  found  in  analogous  fashion  by  differentiating  Equations 
3. 4  and  3.  5  directly  to  get 


(A8) 


(A9) 


Now  uf  must  be  eliminated  and  a  value  wnich  is  a  function  of  the  air  and 
coolant  streamtube  velocity  difference  must  be  substituted  for  e"Xx.  The 
final  velocity  is  known  from  Equation  3.  6,and  a  relation  for  e“Xx  can  be 
derived  by  subtracting  Equation  3.  5  from  Equation  3.  4.  This  gives 


-\x 
=  e 


(A  10) 


Substituting  this  and  the  expression  for  uf  from  Equation  3.  6  into  Equations 
A8  and  A9  gives  the  final  form  of  the  velocity  derivatives  as, 


(All) 


(A12) 
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APPENDIX  B 


1 s _ Outline  of  Non- Equilibrium  Cooling  Computer  Program 

To  analyze  the  cooling  of  a  chemically  reacting  air  stream,  assume 
that  all  chemical  reactions  occur  in  a  streamtube  adjacent  to  a  streamtube 
composed  only  of  SF^  which  remains  a  non-reacting  gas.  Cooling  occurs 
due  to  the  energy  and  momentum  exchange  between  these  two  streamtube  s. 

2. _ Reacting  Streamtube 

The  formulation  for  the  analysis  of  the  reacting  streamtube  is  taken 
from  Reference  5.  In  general,  the  non- equilibrium  state  of  a  gas  is  defined 
by  two  state  variables  such  as  pressure,  P,  and  density,  p,  plus  the 
additional  variables  q^q^,  •  •  •  >*ln  necessary  to  specify  the  internal  non¬ 
equilibrium  of  the  ga&.  The  enthalpy  theij  is  given  by  an  equation  of  state 
of  the  form: 

h  =  h(P,  p,  q,, . . .  ,qn)  (Bl) 

and  the  variables  q^  are  defined  by  a  rate  equation  of  the  general  form 

Dqk 

"EiT  =  wk  P*  ql*  *  *  *  ,qn^  (®^) 

where  Dq^/Dx  is  the  substantial  derivative,  and  the  are  in  general 
complicated  functions  of  the  thermochemical  state.  If  the  rotation  of 
polyatomic  molecules  is  completely  excited  at  its  equilibrium  value  and 
the  effects  of  electronic  excitation  and  ionization  are  negligible,  the 
specify  internal  energy  per  mole  of  a  specie  including  the  effects  of 
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vibration  and  dissociation  i* 


**-v  ♦*».-*». 

1  i  1 


(B3) 


where: 

a 

Cv.  "  sPecific  heat  per  mole  at  constant  volume 

A 

ev.  =  sPecific  vibrational  energy  of  specie 

A 

eD.  =  sPecific  dissociation  energy 
The  total  internal  enerev  of  sn»nV  ; 

.  .  gY  P  ie’  °f  the  reactmg  gas  mixture  is  found 

by  summing  over  all  the  species  to  give 


e  =  Zn>  ( V +  V - 'V )  ■  ™ 

The  temperature  T  is  related  to  the 


(B4) 


of  state 


previous  state  variable  by  an  equation 


p  =  P  L  n.  R  T 


(B5) 


Equation  B5  is  the  same  as  the  familiar  state 
gas  (air) 


equation  for  an  equilibrium 


(B6) 
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with 


*o  = 

By  using  the  definition  of  enthalpy  as 

h  =  e  4-  P/p  (B7) 

Equations  B4  and  B5  can  be  used  with 


(B8) 


which  is  a  special  form  of  the  state  relation  h=h(P,  p,q^}  .  .  .  ,  q^)  applicable 

to  a  mixture  of  reacting  gases.  The  rate  equation  (Equation  B2)  must  now 

be  generalized  into  a  form  applicable  for  the  determination  of  the  species 

n.. 

1 

The  present  computer  program  considers  seven  chemical  species 
02,  N^,  NO,  NO+,  e  ,  O,  N.  The  system  of  ten  reactions  are  represented 
as 
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where  x.  denote  the  chemical  species  listed  in  order  above  and  a  .  and 
1  r,  1 

Pr  •  are  the  stoichiometric  coefficients  of  the  reactants  as  given  in 

Table  B-I.  All  the  reactions  are  listed  in  Table  B-II.  k^  ,k^  are  the 

r  r 

forward  and  backward  rate  constants  which  are  solely  a  function  of 
temperature  in  the  following  way. 


kf  =  A"(T)C  exp^r 


\  =  A‘(T)C  exp4f 


(BlOa) 

(BlOb) 


The  coefficients  of  the  above  two  equations  have  been  taken  from  Reference 
5  and  are  tabulated  in  Table  B —III.  With  this  notation  the  rate  of  change  of 
the  densities  for  the  chemical  species  becomes 


dn. 

i_ 

dissociation 


10 


h  X)  <  @r,  I' "r,  i>  kf  A  lO-i1  r'‘ 
r=l  r  1 


(Blla) 


dn. 

l 


10 


i . 


recombination 


pu  ^r,  i~ar,  i^  ^  (  Pnj) 


r,  l 


(Bllb) 


and  the  net  rate  of  change  of  species  is  then 

d»i  /dn.\  /dn.\ 

dx  y  dx  f  diss.  ~  y  dx  f  recom. 
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The  vibrational  energy  of  the  molecules  is  evaluated  from  the  following 
expressions: 


~  =  g-  P(x)  exp(h.T1//3) 


i=l ,  2 


and  Ty,  the  vibre:  n.onal  temperature  of  the  species,  is  defined  as 


(B14) 


The  constants  0  ,  N  ,  N,  g.,  e.,  h.  for  each  specie  are  shown  in 

i  i 

Table  B-IV. 
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3.  Computer  Program  Inputs  and  Computation  Procedure 


Inputs  to  the  non-equilibrium  flow  cooling  program  are: 


initial  properties  of  air  and  coolant 

V  Tai’  Uai’  Cp  ’  PtOt 

1  1  1  C 


cooling  rate,  and  amount  of  coolant 

^9 


point  at  which  cooling  inititates 


integration  c.ep  «ize 


The  initial  composition  of  the  air  streamtube,  the  total  enthalpy  of  both  the 

air  and  coolant  streamtube s,  and  the  above  derivatives  can  be  computed 

with  these  inputs*  These  derivatives,  together  with  the  derivatives  for 

dn./dx  and  de  /  dx,  allow  the  computation  to  proceed  downstream  in  the 
i 

following  way: 


ha(*i+l>  =  ha<*i>  +  S T  <xi>  Ax 


(B15a) 


aa(xi+l>  =  ua(!ti>+-$r(xi>  Ax 


(B15b) 


u  (x.  .)  =  u  (x.)  + -j —  (x.)  Ax. 
cl  i+l'  cv  i'  dx  1  v  i 


(Bl5c) 
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(B15d) 


de 


'v.(xi+i>  =  ev.<xi> +  ~ar 


dn. 


ni(xi+l )  =  n.(x.)  +  (x.)  Ax 


1  1 


ST 


(B15e) 


Knowing  h  (x)  and  n  (x)  a  new  stagnation  enthalpy  of  the  air  streamtube 


a*  '  a 

can  be  evaluated  as  follows: 


u 


Ha(xHl'=-r(xi+l»  +  ha(xi+l) 


(B16) 


This  also  forms  the  stagnation  enthalpy  of  the  coclant  from  Equation  3.1 

‘  (B17) 


H  (x.^.)  = 
c'  1+1' 


HH  -  H  (x.  ) 

a'  Iti 


Q 


where 

HH  -  (H  +  C~H  ).  .+.  . 

a  Q  c'mitial 

Since  the  velocity  of  the  coolant  streamtube  has  been  calculated,  the 
definition  cf  stagnation  enthalpy  can  be  used  to  find  the  static  enthalpy  of 
the  coolant: 

2 

hc(xi+l>  =  Hc<*i+1>  -  T  (*i+l>  <B18> 

Also,  since  the  coolant  is  assumed  to  be  undissociated  and  to  remain  a 
perfect  gas,  the  temperature  of  the  coolant  is  immediately  found  as 
follows: 

Tc*Xi+P  =  hc^Xi+l^CP  (B19) 
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The  mixing  which  occurs  when  a  coolant  is  injected  into  a  hot  air- 
stream  has  been  postulated  to  take  place  at  constant  pressure.  This 
condition  is,  of  course,  true  at  the  initial  and  final  states  of  the  gases. 
However,  during  the  mixing  process  itself,  the  temperature  of  the  coolant 
and  air  streamtube  are  unequal  ajid  Dalton's  law  of  partial  pressures  does 
not  apply.  Since  the  actual  mixing  process  is  a  highly  complicated  one,  the 
exact  thermodynamic  state  of  the  air  streamtube  is  not  known  and  Dalton's 
law  can  be  postulated  to  hold  as  an  expedient  that  will  allow  the  determination 
of  air  thermodynamic  properties  in  this  region,  knowing  that  the  correct 
final  state  will  be  reached  with  this  assumption.  Thus,  assume  that 


P  =  P 
rOT  air 


+  P 


coolant 


where 


=  p  Z  n.  R  T  + 
*a  1  a 


Pc 


(B20) 


a  =  air  streamtube 

c  =  coolant  streamtube 

m  =  number  of  moles  per  gram  of  species  in  air  streamtube 
R  =  universal  gas  constant 

m  -  molecular  weight 


The  air  and  coolant  streamtube  densities  at  any  point  in  time  is  from 
Reference  1,  as  follows: 

Pc  =  pa  (B21) 
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Upon  substitution  of  the  above,  Equation  B-20  becomes 


[  S  "i  5  Ta  +  1TC-  %  Tc  ]  =  PTOT  <B22> 

From  Equations  B5  andB8,2n.  R  T  can  be  solved  in  terms  of  an  energy 

i  a 

(enthalpy)  relation  and  substituting  in  the  above  equation  gives 


Assume  now  that  mixing  proceeds  at  the 
exchange. 


u 

a 


same  rate  as  the  momentum 


At  the  initial  point  with  no  coolant  injected  into  the  air  streamtube  u  =u  , 

a  3>  j 

and  C.  =  0,  and  the  final  value  of  Cwith  u  =uf  gives 

imtj„l  a  f  6 


C 


f 


Q 

TTC 


Q 


Having  specified  Z  in  terms  of  the  velocity  slowdown  of  the  air  streamtube, 
u 

al 

-  -  1  can  be  substituted  for  C/  1-C  in  Equation  B23,  and  solving  for  the 

a 
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air  streamtube  density 


Once  the  density  and  the  bracketed  terms  are  found,  the  air  stream- 
tube  temperature  can  be  determined  by  using  the  partial  pressure  relation 
Equation  B-20  to  get 


P 

_ a 

Pa  ?  ni  «■ 
1 


(B25) 


and  knowing  T  ,  T  f  u  ,  u  new  derivatives  can  be  formed  for 


dh  du  du  dn.  v. 
a  a  c  1  _ 1 

dx  1  dx  •  dx  9  dx  1  dx 


‘and  the  computation  process  can  be  continued  until  T  =  T  and  u  =  u  at 

a  c  a  c 

which  point  the  reacting  air  streamtube  and  the  coolant  streamtube  will 
have  exchanged  all  the  energy  available  and  computation  stops.  A  typical 
output  of  the  program  is  shown  on  the  pages  along  with  a  listing  of  the 
computer  program  instructions. 
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Non-Equilibrium  Cooling  Computer  Program  Listing 


MAIN-- 


C  MAIN  PROGRAM  FOR  FLOW  OF  A  NITROGEN  OXYGEN  MIXTURE 

COMMON  FMTA,FMTB,FMTC, IDENT ,L INE , X ,H,RO, T, P,U, A, TV.ONOXD. DNOX,DENS 
It ROFREE,PFREE, IM, SMALLR, PO.XO 

DIMENSION  FMTA(96),FMTB(24)» FHTC (12),FMT9(1), THETA VI 2 ) • VNVI2 ) 
DIMENSION  AB{20)»CB(20)(DB(20),AF(20) ,CF(20) ,0F(20) ,CPI20) ,60120) 
DIMENSION  ICB(20),ICF(20)fIA(2if20)fIB(2i,20) ,BMA(20,20 ) , I DENT! 12) 
DIMENSION  DENS(20),EV(20),vm,TV(2),TAUm,EVQ(2ltEXC(2),EXD(2) 
DIMENSION  RKB( 20 1 ,RKF( 20 ) , DNDXDt 20), DNOXRC 20) , DNDXC20) , DEVDXI2) 
EQUIVALENCE  ( (MAD.FIMAD), (IMAX,FIMAX) 

C  SET  UP  CONSTANTS  FOR  CALCULATIONS 

R*8. 31436E7 
THETAVI 1)*2238»97 
THETAV(2)=3336.8 
VNVC 1 )=52. 

VNVl 2 ) “67. 

C  READ  IN  FORMAT  STATEMENTS  FOR  OUTPUT 

READ  INPUT  TAPE  5 , 1001 , ( FMTAI J ) , J=I ,96) , I FMTB(J) , J«I ,24) , (FMTCt J) , 
1  J-l  • 12) 

C  READ  IN  DATA  TO  DEFINE  REACTIONS 

READ  INPUT  TAPE  5 , 1002, I  MAX, IRMAX, ( ABC K) ,CB( K ) , DBt K) ,K* 1, IRMAX) , ( A 
1F(K)»CF(K),DF(K)»K=1» IRMAX) 

C  EVALUATE  COEFFICIENTS  FOR  POWER  OF  T  IN  RATE  CONSTANTS 

00  3  K*i , IRMAX 

ICB!K)=2.0*CB(K)«SIGNF(0.5,CB!K) ) 

ICFtK)*2»0*CFCK)*S IGNF ( 0. 5,CF ( K ) ) 

3  continue 

IFM*IMAX/10 
IFIIFM)  6,6,5 

C  SET  UP  FORMAT  TO  READ  IN  A  AND  B  COEFFICIENTS 

5  I MAD= IMAX*54*IFM 

GO  TO  T 

B6  F IMAD=F I MAX >006000 000000 

B7  FMTD=FIMAD*7400003 10134 

READ  INPUT  TAPE  5»FMTD»t(IA(J,K)»J*l»IMAX),K=l,IRMAX),((IB(J,K)»J« 
11*  IMAX),K*1, IRMAX) 

READ  INPUT  TAPE  5, 1003, ( CPI J ) , J  = I , IMAX ) 

READ  INPUT  TAPE  5, 1004, ( EDI J ) , J«1 , IMAX ) 

C  EVALUATE  COEFFICIENT  FOR  POWER  OF  MASS  OENSITY  IN  RATE  EQUATION 

DO  9  K*i, IRMAX 
I  At IMAX*I,K)=-1 
IB( IMAX*1,K)*-1 
DO  8  J=I , IMAX 
C  CALCULATE  B  -  A 

BMA(J,K)*IB( J,K)-IA(J,K) 

IA( IMAX+i,X) =IAI IMAX*I,K)*IA( J,K) 

IB  t IMAX  +  I,K)  =  IB( I MAX* I ,K ) ♦ IB  I J ,K ) 

8  CONTINUE 

9  CONTINUE 

DO  10  J*3» IMAX 
EVI J ) =0.0 

10  CONI  I NUE 
IM»IMAX-1 

C  READ  IN  JOB  IDENTIFICATION  AND  INITIAL  CONDITIONS 

12  READ  INPUT  TAPE  5, 1005, ( IDENT C J) , J«1 , 12 ) ,X, T,U, P, XC, ALAMOA.CQ, TC,C 
IPC.ANC.HILLNO.tEVl J), J*1 , 2 ) , ( DENS ( J), J*1,IM) 
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Non-Equilibrium  Cooling  Computer  Program  Listing  (Continued) 

MAIN — 


C  INITIALIZE  PRINTEO  AND  PUNCHED  OUTPUT 

WRITE  OUTPUT  TAPE  6,1101, ( I0ENTI J) , J«l» 12) 

LINE-2 

1965  OENSl IMAXI-0.0 

C  CALCULATE  DENSITY  FOR  SUMMATION  NI 

DO  14  J»1,IM 

OENSl  IMAX1-0ENSI IMAX)*DENS( J) 

14  CONTINUE 

RO*P/(T*R«DENS( IMAX) } 

FLCW  =R0*U 
OUMA-O.O 
OUMB-O.O 
CORA— 2.893SE11 

C  CALCULATE  THE  STAGNATION  ENTHALPY 

DO  17  J-l.IM 
DUMA-DUMA+DENSI J) *CP( J ) 

DUMB-DUMB*DENSC J)*(EV( J)-EO( J)» 

17  CONTINUE 

DO  19  J«i,2 

TV(J)*THETAV( J ) /LCGFI R*THETAVl J)/EV< J) *1.0) 

19  CONTINUE 

HT»P*DUMA/(RO*DENS( IMAX ) ) ♦OUHB+O. 5»U*U 
H*HT-0.5*U*U 
UAZERO-U 
80  AHA»1.*CQ 

CPCCOL-CPC*R*ANC 
COOL-C?COOL*CO 
UPRIME  -0. 

NUMBER-0 

C  READ  IN  AN  INCREMENT  CARO 

26  READ  INPUT  TAPE  5, 1006, NSI • IFNI , NI ,OELX 
IF (NS  I )  12,12,1963 

1963  WRITE  OUTPUT  TAPE  6, 1104, NSI , IFNI ,NI «OELX 
LINE  -LINE+2 
IFHEAD»1000*IFNI*NSI-N! 

IF  (IFHEAO)  29,27,27 

C  PRINT  HEADINGS  IF  THERE  WILL  BE  OUTPUT  FOR  THIS  INCREMENT  CARO 

*7  WRITE  OUTPUT  TAPE  6»FMTA 
L INE-L INE*7 
I F t IFNI )  29,29,28 

28  CALL  OUT 

C  SUBDIVIDE  CALCULATION  IF  MORE  THAN  150  STEPS 

29  NLIM-1500-XMODF(1500,NI) 

IFCNSI-NLIM)  30,30,31 

30  NLIM-NSI 

31  NS  I-NSI-NL IM 
NCCUNT-NI 

C  CALCULATE  DERIVATIVES  AND  NEW  VALUES  FOR  ALL  VARIABLES 

00  100  N-1,NLIM 
TCURT-T«*0. 33333333 

TAU(1)-4.924E-5*P*EXPF(0.6464*TCURT) 

TAU(2 )-l. 28E-5-P* EX PF( 0,5493* T CURT ) 

vm-i.o 

Vf  2)- 1,0 
TSO-SORTF(T) 

DC  45  K-l, IRMAX 
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Non-Equilibrium  Cooling  Computer  Program  Listing  (Continued) 
MAIN — 


RKBlK)*ABU)*PONERl  ICBI X ) ,TS0) •EXPFIDBt K) /T ) 
RKF(K)*AF(K) •POWER! ICF ( X ) ,TSQ ) *EXPF (OF I X I /T ) 
DO  43  J*1  *  2 

RXF  I X ) *RKF I X ) • POWER I IAt J,K) ,Vt JJ) 

43  CCNTINUE 
45  CCNTINUE 

00  55  I  *  I »  IM 
00*0.0 
0R*0. 0 

DO  54  X*1 «  IRMAX 
IFIBMA(ItK))  47.54,47 
47  0UMD»BMAt I,K)«RKFIK) 

DUNR-8MAI I,K)*RKBIK) 

DO  51  J*1 , 1MAX 

DUMD*OUMO*POWERI I A ( J, K ) , OENSI J ) ) 
DUHR*DUMR*POWERl IB( J,K) , OENSI J) ) 

51/  CCNTINUE 

LA*IA( 1MAX*1,KJ 
LB*  IB  I  IMAXM.K) 

00*00 ♦OUMO»POWER I  LA ,ROi 
DR*CCUDUMR*POWERlLB,RO) 

54  CCNTINUE 
DNOXD!I)*OD/U 
ONCXR 1 1 ) *DR/U 

ONDXI  I  )*DNOXr.  I 1  )-DNDXRC  I ) 

55  CONTINUE 

00  63  J*l,2 

201  EXA*EXPF( THETAVI J) / T ) 

EVQt J)*R*THETAVI JI/IEXA-l.O) 

DEVDXl J)*(EVQI Jl-EVl J) )«TAU! J) /U 
EVI J) *EVI J) ♦OF VOX  I J)«OELX 

TVC  J)*THETAV(  J)/LOGFIR»THETAVI JJ/EVI  JU1.0I 
63  CCNTINUE 

OENSI  IMAX ) *0. 0 
DUMA-0.0 
0UMB*0.0 
00  65  J*l, IM 

OENSI J)*0ENSIJ)»0N0XIJ)*0ELX 
OENSI IMAX ) *DENS I IMAX ) ♦OENSI J) 

DUMA* DUMA ♦DENS I J) «CP( J) 

DUMB* DUMB ♦DENS I J)*IEVIJ)-E0IJ) ) 

65  CONTINUE 
X*X+DELX 

IFIX-XC)  79,83,83 
83  IF  I NUMBER-l )  82,81,81 
82  HH=HT ♦COOL»TC 
NUMBER*NUM8ER»1 
81  BAH=IU-UPR1ME)/AHA 
UPOX«ALAMBA*BAH 
UPR IME=UPRIME*UPOX*DELX 
DUDX*-UPDX*CQ 
U*U+OUOX*OELX 
HT*H^.5*U«U 
HTPRI M=HH-HT 

HPRIME*IHTPRIM/CQ)-0.5*UPRIMfc*UPRIME 

TPRIME-HPRIME/CPCOOL 
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non 


Non -Equilibrium  Cooling  Computer  Program  Listing  (Concluded) 

MAIN— 


DHDX=-CQ/I l.*CQ*HILLN0)*!CPC00L-U«U/I2.0«AHA*! T-TC II) *ALAMBA« !T-TP 
1RIME) 

H=H*DHDX*DELX 
79  OUTH= ( H-CORR I/7.874E8 

RC=P/ ( (H-DUMB)»DENS( IMAX ) /DUMA*! UA2ER0/U- 1. 0)« ANC*R«TPRIME ) 

PA=P-!UAZER0/U-1.0)*R0*R*ANC«*TPRIME 

T=PA/!RO*R»DENS( IMAX) ) 

A=FLOW/IRO«U) 

NC0UNT=NC0UNT-1 

C  PRINT  OUT  VARIABLES  IF  REQUIRED 

IF1NC0UNT)  71,71,100 
71  CALL  OUT 
NCCUNT=N1 
100  CONTINUE 

IF( IFNI*(NI-NCOUNT) )  105,110,105 
105  CALL  OUT 

110  IFUFHEAD)  120,115,115 

IF  THERE  HAS  BEEN  OUTPUT  FOR  THIS  INCREMENT  CARD  PRINT  OUT 
VALUES  OF  VARIABLES  FOR  RESTART  OF  CALCULATIONS  AND  INITIALIZE 
PRINTED  OUTPUT 

115  WRITE  OUTPUT  TAPE  6, 1 102, ( E V( J » , 1, 2  I , FLOW.OUTH, TPR IME.UPRlME, HPR 
1 IME 

WRITE  OUTPUT  TAPE  6, 1 101 , ( IDENTI J > , J=1 , 12 ) 

LIKE*2 
120  IF(NSI)  29,26,29 

1001  FORMAT! 12A6) 

1002  F0RMAT(2I5/CE20.5,F10.1,F20.0> ) 

1003  FORMAT! 10F7. 2) 

1004  FORMAT ( 6E 12. 3) 

1005  FORMAT! 12A6/4E15. 5/3E15. 5/4E15.5/2E15. 5/ !4E15. 5 1 ) 

1006  FORMAT! 315, 1PE15. 5) 

1101  FORMATC43H1N2TROGEN-OXYGEN  FLOW  THROUGH  A  STREAM  TU8E/1X, 12A6I 

1102  FORMAT ! 2 7H0V IRRATIONAL  ENERGY  OF  C2  =, 1PE 15. 7/27H  VIBRATIONAL  ENER 
1GY  OF  N2  *, 1PE15.7/17H  MASS  FLOW  RATE  =  , 1  PEI  5. 7/ 1 1H  ENTHALPY  «,1PE 
215.5/9H0TPRIME  =,1PE15.5/9H  UPRIME  »,1PE15.5/9H  HPR I ME  *,1PE15.51 

1103  FGRMAT(llA6,6X,A6,lH0,m 

1104  FORMAT! 1 5H0I NCREMENT  CARD, 31 5 , 1PE 10. 2 ) 

ENDI1, 1,0, 0,0, 1,1, 1,0, 0,0, 0,0, 0,0) 
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TABLE  B-I 


STOICHIOMETRIC  COEFFICIENTS 


d 
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N2 
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1 

0 

0 

3 

1 

0 

0 

0 

0 

0 
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1 

4 
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1 
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M  =  2  n. 
i 
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reactions 
table  b-ii 


d°2 

dx 


"  pu  j^kf^P0)  "^(pO)  (p02)-kf^(p02)(p0)+kb^{p0)3-kf  (pOz)(pM) 
+kb3(P°)2(PM)-kf^(p02)(pN)+kb^(pN0)(p0)-k^(pN2)(p02)+kb  (pNO)2J 


dO 

3F 


(P02)2-2kbi{p0)2(p02)+2kf2(p02)(p0)-2kb2(p0)3+2kf3(p02)(pM) 


dN-, 

Ld 

■sr 


=-l  K 
pn  L  fi 

-2kb3(p°)  (pM)+kf^(pNO)(pM)-kb^(pN)(pO)(pM)+kf^(pN)(p02)-kb  (pNO)(pO) 
-k{  (pOHPNz)+kh  (pNO)(pN)-k  (pN>(pO)+k,  (pNO+)(pe-)l 

0  8  10  Dio  /  J 

"  pu  ['kf9(pN2)  +kb9(pN2)(pN)2_kf5<PN2HpM)+kb^(pN)2(pM)  [ 

” kf 9 ( PN2 } ( P°2 ) +kb  (PN°)2] 


dN 


Pu  [+2kf4(pN2J  -2kb4(pN)2{pN2)+2kf  (pN2)(pM)-2kb  (pN)2(pM) 

5  5 

+kf^(pNO)(pM)-kb^(pN)(pO)(pM)-kf^(pN)(pOz)+kb  (pNO)(pO) 
+kf8(pO)(pN2)-l^8(pNO)(pN)-kf^(pN)(pO)+kb^(pNO+)(pe.-)j 


dNO  _  de 
~3x 


dNO 

~dx 


ir  |kf  (PN)(P°)"kb  (pNO+)(pe')l 

L  10  Dio  J 

"pH  ["kf6(pN0^PM)+kb6<PN)(P°)(pM)+kf^(pN)(p02)-kb  (pNO){pO) 
+kf8(P°)(PN2)-kbg(pNO)(pN)+2kf^(pN2)(p02)-2kb  (pNO)2j 
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